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Abstract: We analyze a class of bottom-up holographic models for low energy thermo-electric transport. 
The models we focus on belong to a family of Einstein-Maxwell-dilaton theories parameterized by two 
scalar functions, characterizing the dilaton self-interaction and the gauge coupling function. We impose 
spatially inhomogeneous lattice boundary conditions for the dilaton on the AdS boundary and study the 
resulting phase structure attained at low energies. We find that as we dial the scalar functions at our 
disposal (changing thus the theory under consideration), we obtain either (i) coherent metallic, or (ii) 
insulating, or (iii) incoherent metallic phases. We chart out the domain where the incoherent metals 
appear in a restricted parameter space of theories. We also analyze the optical conductivity, noting that 
non-trivial scaling behaviour at intermediate frequencies appears to only be possible for very narrow regions 
of parameter space. 
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1 Introduction and Outline 

The holographic AdS/CFT correspondence provides a valuable theoretical framework for exploring the 
dynamics of strongly coupled field theories. In recent years this approach has been found to be quite useful 
in building qualitative models for a variety of phenomena encountered in many-body physics, ranging 
from superconductors and superfluids, to non-Fermi liquids [18, 26, 29]. The general thrust of the initial 
discussion has been to consider relativistic quantum field theories which are deformed by relevant operators 
to access non-trivial phases in the infra-red (IR). In interest of computational simplicity, much of the 
discussion was focused on situations where spatial homogeneity is retained. 

Given the progress in these homogeneous models it is natural to consider less symmetric systems as 
a means to accessing more realistic scenarios. For instance, consider the linear response functions which 
may be efficiently computed in the holographic context. In the absence of dissipation the low frequency 
behaviour of correlation functions often display IR divergences. While some of these can be cured by 
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examining the theory at hnite temperature (effectively using a thermal IR regulator), there are some issues 
that are incurable by this simple approach. 

A simple case in point is the physics of electrical conductivity. Naively, we might try to compute 
this quantity in a translationally invariant system by monitoring the response to turning on an external, 
time dependent electric field. Using the Kronig-Kramers relations one would then find that momentum 
conservation results in a divergent zero frequency value of the imaginary part of the conductivity. As 
a consequence recent efforts have been focused on constructing models with more realistic conductivity 
behaviour. This in turn is motivated by the need to better understand the experimental observation of 
metal/incoherent-metal/insulator transitions in the phase diagram of high temperature superconductor 
materials [4]. 1 

The first models for realistic behaviour of (thermo)-electrical conductivities in holography was accom¬ 
plished in [27, 28]. These authors considered situations where the spatial homogeneity is explicitly broken 
by a background lattice of sources, and studied thermo-electric transport in the resulting holographic dual 
background. By use of the lattice, they were able to demonstrate that the low frequency conductivity 
was hnite, and well fitted in their set-up by the Drude form. More curiously, their analysis revealed a 
mid-range scaling of the AC conductivity with frequency, which furthermore agreed with experimental 
results in cuprate systems [34], 

Since this seminal work, various groups have attempted to understand the thermo-electric properties 
of holographic systems. One of the key aspects to understanding this phenomenon is to ascertain the 
potential low energy behaviour of strongly coupled systems subjected to inhomogeneous sources for relevant 
operators. A classification of emergent IR phases crucially provides one with a picture of what to expect 
for the transport properties of the system at low energies. Extensive work has already been carried out in 
characterizing the possible IR geometries which may emerge when translational invariance is relaxed; see 
for example [10, 15] where the IR flow of perturbative inhomogeneous modes was analyzed. We now have 
a reasonable understanding for the spectrum of possibilities that can occur in holographic systems, with 
low energy behaviour ranging from metallic (with various dressing) to insulating. Part of the motivation 
for the present work was to get a better picture of the landscape of possibilities within the framework of 
bottom-up holography. 

We should also note other innovative approaches towards modelling aspects of inhomogeneous systems. 
These include the techniques of massive gravity [1, 5, 6, 9, 35] and Q-lattices [2, 11-13, 16, 17] (cf., also 
the memory function approach [31, 32]). The Q-lattices are particularly appealing as, though translational 
invariance is broken, the metric remains homogeneous and the equations of motion remain relatively 
tractable ODEs. Unfortunately the fine tuned nature of this symmetry raises the possibility that some 
of the behaviour observed in the resulting solutions may be equally non-generic. Likewise while much 
progress has been made using the massive gravity approach, questions remain about the interpretation of 
such models in the dual field theory. 

With these issues in mind let us turn back to explicit breaking of translational invariance in one of the 
spatial directions of the dual held theory, as described in [27, 28] (see also [30]). As mentioned above, the 
introduction of the holographic lattice and solutions of the resulting PDEs mitigates the zero frequency 

1 A nice summary of the issues involved in characterizing such behaviour in many-body systems can be found in [20]. 
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delta function in the conductivity, leading to a conventional Drude form for low frequencies. An extremely 
thorough examination of the subject in [14] combines technology gained from Q-lattice calculations with 
the numerical approach of [27, 28, 30] to examine the electric, thermoelectric and thermal conductivity. 
Among the results of this analysis were the confirmation of the Drude regime for the AC conductivity and 
the existence intermediate resonance peaks. It was also noted that in the model under consideration no 
evidence of a rnid-IR scaling regime was discovered. 

In this note we aim to build on the work of [14, 28, 30] by exploring a particular class of Einstein- 
Maxwell-Dilaton models. These models are characterized by two functions of the scalar dilaton field <f> - 
these parameterize the scalar self-interaction through a scalar potential and a gauge coupling function, V (4>) 
and Z(4>) respectively. This parameterization allows us to access a range of effective field theories which 
are of interest in various condensed matter applications. For example, examination of the homogeneous 
solutions has shown that these models are good candidates for metal/insulator transitions upon breaking 
of translational invariance [8]. The range of allowed IR phases can be understood in terms of an effective 
near horizon potential, obtained by approximating E(3>). One finds that this potential modulates the 
spectrum of ingoing excitations near the horizon from continuous (metallic) to discrete (insulating) as a 
function of parameters. Qualitatively, the physical features of these models can be understood by noting 
that scalar potential interpolates between zero at the conformal boundary and a runaway behaviour in the 
near-horizon region. Similarly the gauge coupling starts at some finite value near the conformal boundary 
and undergoes exponential growth or decay towards the horizon. 

We will examine the electric conductivity of these models as a function of parameters. We employ two 
lines of attack, both of which involve first numerically finding the bulk solution to sourced inhomogeneities 
in the UV. This holographic lattice solution can then be explored for transport. The simplest analysis is 
then to extract the DC conductivities via a membrane paradigme-esque formula using the techniques of [14] 
(cf., §3). This approach has the advantage that we only need knowledge of the bulk solution and therefore 
gives us a simple way to probe the IR phase by analyzing the DC conductivity [28]. Computing the 
AC conductivity however requires that we also have the solution of the linearized perturbation equations 
around the numerically constructed background [28] . We carry out this exercise to obtain the full frequency 
dependent transport, in the process using the aforementioned DC conductivity as a non-trivial check on 
the results we obtain. We discover the existence of metal-insulator transitions as a function of parameters 
and also, in finely tuned cases, the potential presence of a rnid-IR scaling regime. We then examine the 
persistence of these behaviours as a function of parameters. We find that that the phase changes we 
discover via monitoring the behaviour of the response functions correspond to inhomogeneity mediated 
changes in the form of the scalar potential in the IR. 

The outline of the paper is as follows. In §2 we present a basic review of the holographic set-up, 
pausing to note the ingredients we pick in our model and their potential effects on the low energy dynamics 
of the system. We also give a short synopsis of the numerical scheme we employ to study the physical 
transport. In §3 we revisit the arguments of [14] to directly extract the zero frequency conductivities in 
terms of horizon data. In course of this analysis we take the opportunity to explain the relation between 
the various conductivities in the hydrodynamic limit (which corresponds to w < T and a suitably dilute 
lattice). We illustrate that in this regime there is a single transport coefficient which can be taken to be the 
electrical conductivity. In §4 we present our results for the various transport coefficients. We also comment 
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on the various phases and scaling regimes which we observe as we explore a representative subspace of 
configurations for the scalar potential and Maxwell coupling. We conclude with some discussion in §5. 
Some technical points about our numerics are collated in Appendix A. 


2 Setup 

We start by outlining the holographic set-up we use to explore the physics of conductivity in 2+1 dimen¬ 
sional quantum critical systems. We explain the basic ingredients we employ in our phenomenological 
modelling as well as the salient features of the control functions we introduce in parameterizing the holo¬ 
graphic system. Following this discussion we go on to describe the numerical scheme used to construct the 
gravitational solutions of interest and the computation of the physical conductivities therein. 


2.1 Background 

We take a bottom up, phenomenological approach to holography, and following [14, 28, 30] we choose a 
model with the minimal ingredients necessary to calculate conductivities in a 2+1 dimensional quantum 
critical theory. We assume that the field theory is holographically dual to gravitational dynamics in an 
asymptotically AdS 4 spacetime. The model for the holographic dynamics is simply Einstein-Maxwell 
theory, which we couple to an additional neutral scalar field, with potential term, V(3>), and a gauge 
coupling function, Z(4>). The scalar field allows us to discuss models with more general behaviour in the 
IR than the local criticality characterized by the Reissner-Nordstrom-AdS 4 (RN-AdS 4 ) black hole. This 
results in an action of the form (setting ^AdS = 1 ) : 

5 = Jd 4 xV^g (R + 6-\z($)F ab F ab -V a <S>V a <S>-^V(<S>)^ . ( 2 . 1 ) 

The dynamical equations of motion which we will solve are then simply 


V a (Z(<S>)F b a ) = 0, 

V a V a <b — V'(<1>) = 0, 

9abF cd F cd ^j 


Rab + 3 g a b ~ ( + - 


2 g ab v(<s>) + z(<s>) 


/? /?+ _ 

r ac r b ^ 


= 0 . 


( 2 . 2 ) 


The field theory will be taken to live in Minkowski spacetime and we shall thus work with the conventional 
parameterization of AdS geometry in Poincare-like coordinates. The radial coordinate is taken to be z € M+ 
and the conformal boundary located at z = 0. 

The potential and the gauge coupling functions are our control functions which allow us to modulate 
the IR dynamics. In [ 8 ] these were parameterized to be of the form F(4>) = W' 5 '*’, and Z(<&) = e 7 ^ 
respectively. It was argued that these choices allow for a range of locally critical IR behaviours as one 
scans over the parameters { 7 , 5}. As we are interested not only in the IR behaviour, but also in translating 
the local critical dynamics therein onto the AdS boundary, we need to ensure that any such IR geometry can 
be patched to the asymptotically AdS 4 region. We must therefore generalize the form of these functions. 
In doing so we choose to fix the potential to have a Taylor expansion around the origin of field space of 
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the form V($) = <h 2 + • • •. This choice corresponds to an effective conformal mass term m 2 = —2 for 
the scalar ensuring that we have simple fall-offs (with non-normalizable and normalizable being 2 and z 2 
respectively asymptotically). A convenient choice for the functions which respects these constraints turns 
out to be 

V($) = 4(1 - cosh(u $)) (2.3) 

yZ 

Z($) = e v<s> 

where we have chosen to focus on a one-parameter family of theories, parametrized by v. This corresponds 
to the case 7 + 5 = 0 in the notation of [ 8 ]. We leave exploration of the full parameter space for future 
study, and focus on this subset henceforth. 

For computational simplicity we introduce inhomogeneity in our model by choosing to break transla¬ 
tional invariance in one-direction of M 2,1 as in [14, 28, 30]. We construct a lattice in the x-direction while 
maintaining translational invariance in the y-direction . 2 The x-translation breaking boundary conditions 
are imposed by choosing an inhomogeneous normalizable mode for the scalar field of the form 


<Fi(x) = C cos (kx ), 

where C is the amplitude of the inhomogeneity and the k the wavenumber of the lattice. 

The manner in which this sourced inhomogeneity deforms the near horizon geometry from its homoge¬ 
neous behaviour, and the effect that this has on the conductivity is the principal focus of this work. From 
the dual field theory point of view this corresponds to determining how the UV parameters of the theory 
change the trajectory of the renormalization group flow to create different phases of matter in the IR. As 
we explain below the presence of radially conserved quantities in the theory mean that a great deal about 
the linear response functions, and therefore the phase of the dual field theory, may be extracted simply 
from knowledge of the near horizon geometry. 

Given these boundary conditions, a suitable metric ansatz for the investigation of these inhomogeneous 
phases is [28]: 

ds 2 = 4 Qtt dt 2 + dz 2 + Q xx (dx + z 2 Q xz dz) 2 + Q yy dy 2 ^ (2.4) 

with the functions Q a b(z,x ) depending both on z and x, thanks to the inhomogeneity. In addition we 
introduce a redshift factor f(z) which, upon exploiting the scaling symmetries in the problem to fix the 
horizon to be at z = 1 , can be chosen to be: 


f(z) = (1 - z)P(z) = (1 - z)( 1 + z + z 2 — Mitt) 


(2.5) 


The factor can be thought of as a convenient parameterization of the temperature, T ; they are related 
via 


T = 


P( 1) 6-y 2 


47T 


8tt 


( 2 . 6 ) 


" We use the word “lattice” loosely for we choose not to impose any commensurability conditions between the charge 
density and the unit cell. 
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This choice is useful as in the homogeneous limit the standard RN-AdS 4 black hole may be recovered by 
setting Q t t = Qzz = Qxx = Qyy = 1 , Qxz = $ = 0 , A = (1 — z) g, and g = g\. We choose to measure all 
physical quantities relative to a fixed relation between T and g which we fix by taking g = g 1 throughout. 
Also, for convenience, we choose A = (1 — z)Aq(z,x) dt ensuring thereby that the timelike component of 
the gauge field vanishes at the horizon. 

In order to ensure that our problem is well posed we must supply both consistent boundary conditions 
and an appropriate gauge condition to remove the gauge redundancy of the Einstein equations. We choose 
to work in DeTurck (or harmonic) gauge - this is achieved by modifying the Einstein tensor via the addition 
of a new term involving the so-called DeTurck vector field £ a : 

G5 = G ab - V (o &, (2.7) 

«“ = s“‘ra»)-rs,(s)) 

Here G ab is our original Einstein tensor, £ a is the DeTurck vector, and G^ b is the modified tensor appearing 
in the DeTurck equations. The DeTurck vector is defined using the difference in the Christoffel symbols, 
r, associated with our metric of interest, g, and a suitably chosen reference metric, g. The reference metric 
should have the same asymptotic and conformal structure as the metric we are attempting to solve for. In 
our case we have found it convenient to use the RN-AdS 4 metric as the reference metric. It can be shown 
that for the metric ansatz (2.4), the DeTurck equations are elliptic [25] and therefore can be solved as a 
boundary value problem. 

As we are interested in solving the original Einstein equations, we must ensure that the DeTurck vector 
vanishes on-shell. Thus we choose boundary conditions such that the the DeTurck vector vanishes on the 
boundary. Provided the solutions to our problem are unique and smoothly dependent on the choice of 
boundary conditions, this should ensure the vanishing of the DeTurck vector. In practice, we check that 
the DeTurck vector is zero to a high numerical precision, so we may be confident in the veracity of our 
solutions. 

After the gauge fixing procedure has been completed we obtain seven independent equations in the 
seven unknowns, (Qtt, Qzz, Qxx, Qyy: Qxz, A)). It can be checked that taking appropriate linear combi¬ 
nations of the equations decouples their principal parts, and each has an elliptic form. 

We next turn to the discussion of the boundary conditions. At each boundary we require one boundary 
condition for each of the seven dynamic fields. We choose the following conditions at the conformal 
boundary: 


Qtt( 0, x) = Q zz ( 0 , x) = Q xx ( 0 , x) = Q yy ( 0 , x) = 1 , Q xz ( 0 , x) = 0 ( 2 . 8 ) 

4>'(0, x) = d-i(x), A 0 ( 0 , x) = g 

The motivation for the metric boundary conditions is simply that one obtain AdS 4 spacetime at z = 0. 
Labeling the Dirichlet boundary condition on Ao(0,x) as g is consistent with its interpretation as the 
chemical potential in the dual field theory. The condition on the scalar field fixes the non-normalizable 
mode to be inhomogeneous and thus sources inhomogeneity in the system. 

In the IR we must first impose the regularity of the black hole horizon. This is done by requiring 
Qtt(l,x) = Q zz (l,x ) which ensures the surface gravity is constant along the horizon. The remaining 
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boundary conditions may then be found by expanding the equations of motion to leading order in (1 — z). 
As two of the equation expansions are degenerate at this order we obtain the correct number of horizon 
boundary conditions. An intuitive understanding of these conditions can be found by solving these leading 
order expansions for the radial derivatives of the fields at the horizon. We then see that these equations fix 
the radial derivatives of the fields in terms of their horizon values, and first and second spatial derivatives 
of the same fields. Therefore our horizon conditions consist of one Dirichlet condition, necessary to define 
a regular horizon, and six (mixed) Robin-type conditions. 

At the spatial directions we impose that all fields are periodic, with period L = This allows for 
the sourced inhomogeneity to be a linear combination of the basic harmonic cos (kx), and any of the higher 
harmonics cos (nkx) for any integer n. For example, [14] and, in a different context, [3, 23] have constructed 
solutions with sourced inhomogeneity consisting of multiple modes with random relative phases in order to 
represent a ‘dirty’ lattice. While such investigations are interesting we postpone the construction of such 
solutions and the investigation of their thermodynamic and transport properties to future work, and focus 
below on constructing solutions sourced by the single harmonic cos (kx). 

2.2 Perturbations and linear response 

We now turn to analysis of linearized perturbations, needed to extract information about the conductivity 
of the QFT in the linear response regime. We will primarily be interested in computing the AC (optical) 
electrical conductivity as a function of frequency and temperature. The results have an intrinsic interest, 
and can also provide a point of comparison to the DC conductivity formula which we will discuss in §3. 
They thus provide a non-trivial check on our results. Combined with the DC conductivities, this provides 
a comprehensive picture of the behaviour of the theory at the temperature regimes we probe. 

The calculation of the conductivity in a holographic theory entails solving the linearized perturbation 
equations derived from Eqs. (2.2). To derive these we perform the following expansion of the fields: 

9ab — 9ab T e h a b, A a — A a T e b a , 4? — T c tj (2.9) 

with h a b , b a and rj being the metric, gauge and scalar perturbations, respectively, and the hats indicating 
background fields. We work in leading order in e and use as the background the solutions to Eqs. (2.2), 
subject to the boundary conditions described in §2.1. The symmetries of our background allow us to 
set by,hty,h Z y,h xy to zero. Therefore the non-trivial components our metric and gauge perturbations are 
restricted to {htt,h zz ,h xx ,h yy ,ht x ,ht z ,h zx ,b x ,b y ,b z ). In addition since our background is static we may 
Fourier decompose the time dependence of the perturbations: 

h a b = hab{z,x)e~ lut , b a = b a {z,x)e~ zujt , T] = ij(z,x) e~ lU}t (2-10) 

As our equations are linear the time dependence is encapsulated in the factors of the frequency u present 
in the equations. 

In analogy to the background equations we must impose appropriate gauge constraints on our problem 
in order to remove spurious degrees of freedom and obtain a match between the number of fields and 
independent equations of motion. In the case of the perturbation equations the gauge redundancy is more 
complex as, in addition to the diffeomorphism invariance of the gravity sector, the U( 1) invariance of the 
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gauge sector must also be taken into account. In order to remove these gauge redundancies we choose to 
impose the deDonder and Lorentz gauge conditions: 

n = V“ ( h ab - h 9 ^) =0, x = V a 6 a = 0 (2.11) 

The procedure we use to impose the gauge conditions is very similar to that used to form the DeTurck 
equations in equation (2.7). For example, in order to impose the deDonder gauge we use a gauge fixing 
term constructed from a gauge transformation, r b . Using the expression for the variation of the Ricci tensor 
and subtracting the gauge fixing term V= V( a V c /i^ — |V a V;,/i gives us: 

SRab - V (a T fc) = - Vv c h afe + \v c V a ht + \VcVbK - \v h V c h c a - ^V a V c hg (2.12) 
= -\v c v c h ab + i[V c , V a ]h c b + ^[V c , V b }h c a 

Thus we see that the principal part of the equation is hyperbolic and the goal of the gauge fixing procedure 
is achieved. Similarly, in the case of the gauge field equations we add a gauge fixing term of the form 
a(z, x)V b x- Here a(z,x ) is some combination of the background fields which is determined by requiring 
that the principal part of each of the gauge equations takes the appropriate hyperbolic form. 

The end result of the gauge fixing procedure is 11 independent equations in 11 dynamic fields whose 

principal parts are decoupled and of a hyperbolic form. We can now consider the boundary conditions. At 

the conformal boundary, these are simple - we require b x to source the external external electric field and 
all other perturbations to vanish suitably quickly such that asymptotically AdS .4 is maintained. Therefore 
we choose the following: 

h tt { 0,x) = 0, h zz ( 0,z) = 0, h xx (0,x) = 0, h yy (0,x) = 0, h tx (0,x) = 0 (2.13) 

h tz { 0 ,x) = 0 , h zx ( 0 ,x) = 0 , b t (0, x) = 0 , b z (0 ,x)=0, b x (0,x) = l 


Here we have used the linearity of the equations to choose a convenient scale for the magnitude of the 
external electric field. 

The situation at the horizon is more complicated. Physically our boundary conditions must reflect the 
fact that all excitations are in-falling at the horizon. This is most readily observed by passing over to (the 
regular) ingoing coordinates Eddington-Finkelstein coordinates and ensuring that both the stress-energy 
and Einstein tensors are regular at the horizon. This will determine the leading scalings of the fields at 
the horizon. Our results, which of course agree with those of [28], are: 


h tt (z,x ) = V{z) h r t ^ 9 (z,x), 
h xx {z,x) = V(z) h r xx 9 (z,x ), 

htz{z,x) = h r t l 9 {z,x), 

h Z z(z,x ) = h r z e z 9 (z,x), 

b t (z,x) = V{z) b 7 t e9 {z,x), 
fj(z, x) = V(z) fj re9 (z , x) 


hyy(z,x) = V(z) h r yy 9 {z,x ) 

h tx {z,x ) = V{z) h r tx 9 {z, x) 
h xz (z, x) = h r xz 9 {z, x) 

b z (z,x)='^ z ~b r z e9 (z,x) 
b x (z,x) = V{z) b r x e9 (z,x ) 


(2.14) 






where 


V(z) = (l-z)~%, tr> = £ (2.15) 

captures the leading non-analytic behaviour of the fields near horizon and the remaining analytic part is 
indicated by the superscript “reg”. Those regular fields may be expanded in a power series in (1 — z) 
near the horizon. Including the leading and subleading orders in this expansion, one identifies 15 possible 
non-degenerate boundary conditions - 4 from the leading order and 11 from the subleading order. Only 11 
of these, chosen for numerical convenience, are imposed. Consistency of the equations of motion demands 
that the remaining 4 conditions vanish on-shell. The vanishing of these constraints provides a non-trivial 
check on the accuracy of our numerical solutions. Finally, we will also demand that the solutions obey 
appropriate periodicity conditions in the x-direction. 

Once the solutions of the perturbation equations are available the AC conductivity can be calculated 
as [14, 28] as: 

o{uj, x) = X , with b x (z, x) = 1 + j x (x)z + 0(z 2 ) (2-16) 

i uj 

From this the DC conductivity may be extracted by taking the uj —» 0 limit of the Re(u). We now turn to 
describing an alternative way of extracting the DC conductivities, involving knowledge of the background 
solution alone. 


3 Analytic Expressions for the DC Conductivities 


We now describe how to extract information regarding the electric, thermoelectric, and thermal conduc¬ 
tivities from knowledge of the background fields alone. This discussion closely follows that of [14] (adapted 
to our coordinates). 

We begin by briefly reminding the reader of the form of the linear response transport equations in a 
2 + 1 dimensional theory. At finite chemical potential, and therefore finite density, the heat and electric 
currents may mix. Therefore Ohm’s Law takes the more general form of: 



/ a aT 
\aT R 



(3.1) 


In our case J = J x is the electric current and Q = T tx — / iJ x is the heat current and T = — \7 X T 
is the thermal gradient. Our approach will be to perturb our holographic system such that E and T 
are introduced, consecutively, in the dual field theory. We will then find that certain radially conserved 
currents will allow us to read off the conductivity pairs of (a, a) and ( R , a) from the near horizon geometry 
of the background. 


3.1 Response from hydrodynamic perspective 

At the outset however, we should remark the following. The DC response of the system comprises of two 
components: a response due to impurity scattering which leads to momentum relaxation and Drude peak 
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behaviour and a more primitive hydrodynamic response present even in translationally invariant systems. 
Specifically, we may write as in [24] the low frequency conductivity in the form 

= , KaT + &h(T, p ), (3.2) 

1 — 1(jJT 

to emphasize that there is a calculable contribution to transport even when momentum dissipation is 
swtiched off. Of course, this contribution has to be extracted after subtracting off the divergent DC 
conductivity arising from the delta function contribution at u = 0 in the r —> oo limit. 

In the hydrodynamic limit 3 there is only a single response coefficient given by ah- To see this, note 
that the limit involves frequencies and momenta which are much smaller than the characteristic thermal 
scale. The lattice if present is treated as a spatial long-wavelength perturbation about a homogeneous 
background. The hydrodynamic energy-momentum and charge currents to first order in spatio-temporal 
gradients take the form 


T^ u = eu^u" +pP^ - 2 n(T, p) a ^ - C (T, p) 0 P^ 

J» = pu» + a h (T,p) (3.3) 


where P ab = g ,iu + u v , is the spatial projector, a^ w is the shear tensor and 0 is the fluid expansion. 
The thermodynamic data is encoded in the energy density e, pressure p and charge density p and we have 
assumed that the underlying system is relativistic with Lorentz invariance being only broken by the choice 
of inertial frame picked by the fluid (through rffi). 4 

The shear and expansion contributions are tensor and scalar modes in the hydrodynamic expansion, 
leaving the contribution coming from the electric field and the gradient of the chemical potential and 
temperature to be the only vectorial part of transport. The existence of a single vectorial transport 
encoded in ah is related to the fact that fluids are required to satisfy the second law of thermodynamics. 
Allowing for three independent vector transport coefficients is inconsistent with the existence of an entropy 
current with non-negative definite divergence, cf., [24], Translating this observation we require that in the 
hydrodynamic limit: 5 


& — &h i 


T 




k = —a p = 



(3.4) 


This can be used to check some aspects of numerics in the high temperature regime (where all our models 
exhibit metallic behaviour). 

3 We use the phrase ‘hydrodynamic limit’ to refer to the low energy description of translationally invariant systems which 
is traditionally well described by relativistic hydrodynamics. Our aim here is to illustrate the fact that underneath the Drude 
peak, there is a single frequency independent ‘hydrodynamic conductivity’ captured by ah- 

4 In writing this expression we have chosen to fix some field redefinition ambiguity inherent in hydrodynamics by demanding 
that be a timelike eigenvector of the stress tensor with eigenvalue being the energy density (Landau frame). 

5 A version of this relation was derived without invoking the second law but using holography in [21]. They however derive 
a relation that mixes zeroth and first order in gradients, which is at odds with interpreting hydrodynamics as a low energy 
effective theory. Disentangling this leads to the relation we quote in (3.4). 
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3.2 Membrane paradigm for response 

Having understood what the relations we expect are, we can now turn to asking what the holographic 
modeling has to say about the transport coefficients of interest. 

We start by considering the electric and thermoelectric conductivities, a, and a. To extract the 
DC conductivities we modify the perturbation ansatz: (hft, h zz , h xx , h yy . ht x , hf Z . h zx , b y ,b z , rj) are time- 
independent and b x = b x (z,x ) — eEt . We may then proceed as described in §2.2 and derive the corre¬ 
sponding equations of motion. 6 It may then be easily checked that the linearized perturbation equations 
for the gauge field imply d x (y/^g Z F xr ) = 0 and d r ( y / zr g Z F rx ) = 0, and therefore' 

J = \J—g Z F xr = constant. (3-5) 


As this quantity is a constant it can be evaluated anywhere including at the horizon. In order to derive 
an explicit form for such an expression we wish to extract the leading scaling of the fields near the horizon. 
This may be done, as above, by transforming into ingoing coordinates and demanding the regularity of the 
stress-energy and Einstein tensor. The results are as follows: 


h t t(z, x) ~ 47r/i 6 t (x)(l - z) + 0(1 - z) 2 , h tz (z,x) ~ h° tz (x) + 0(z - 1) 


Qxx (1)#) 


h tx (z,x) ~ Y + °( 1 - h zx(z,x) ~ h rx( x ) + °( 1 - z ) 

h xx (z, x) ~ h xx (x) + 0(1 - z), hyy (z i x) ~ hy y (x) + 0(1 - z) 
rj(z,x ) ~ rf(x) + 0(1 — z), bt(z,x) ~ b®(x) + 0(1 — z) 


b z (z,x ) ~ 


1 


47tT(1 — z) 
and with the restrictions that: 


b° z (x) + 0(1 - z), b x (z,x) ~ log (47rT(l - z))b x (x) + 0(1 - z) 


E 


h° zx (x) = h° tx (x ), h° zz (x) = 2 h° tz (x) - h° tz (x ), 6° (x) = —, b° z (x) = 6? (x) 


(3.6) 


(3.7) 


in order that Einstein and stress energy tensor be regular in this coordinate system. Note that for conve¬ 
nience we have chosen to include factors of the background metric fields evaluated at the horizon in our 
definition of h ^ x , h® x . 

Expanding our expression for J and utilizing these near horizon expansions we obtain the following 
equation: 



- A 0 



+ b° t '(x) - E = 0 


(3.8) 


where all background fields have been evaluated at the horizon. 

6 For the purpose of the argument that follows the details of gauge-fixing of the equations is inessential. 

' We will use J and Q to denote the bulk conserved quantities, which will of course agree with the boundary charge and 
thermal currents. 
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The next step is derive a similar expression for another conserved quantity which, as explained in [14], 
corresponds to the thermoelectric conductivity. To do so we note that if £ is a Killing vector satisfying 
C^F = 0 then we may define: 

gap = + 1 £[app]a + I (^ _ 2 0) F ap (3.9) 

with ip and 6 defined as C^A = dip and i^F = dO}' This has the important property that 

V a G ap = 3^ (3.10) 


Choosing the Killing vector, £ = dt and utilizing the above identity and the equations of motion it can be 
shown that V x Q xr = d x (^F r g Q xr ) = 0 and '\7 r Q rx = d r {y/—g Q rx ) = 0 and that therefore that 8 9 

Q = \J—gQ rx = constant. (3-11) 


If we now choose that 6 = — (— eEx + Aq(z, x) + e bt(z, x)) and ip = eEx, and expand Q around the 
horizon as we did J , we obtain at leading order: 

Q = — 4-7T T ht x (x) . (3.12) 


This indicates that h® x (x) = h° tx is a constant. Working at subleading order we obtain a constraint 
expression which, with a little rearranging, can be written as: 



h° tx 

2-jtTQzz 


nTdz log 


Qxx Qzz 

Qtt Qyy 


8ttT + 6 


J ^0 A , 0 h l{x) 
4vr T Q zz ) dx Q zz 


(3.13) 


where again all background fields are being evaluated at the horizon. Integrating this expression over one 
period of the background, the last term vanishes as a total derivative whose boundary contributions cancel 
as a result of periodicity. We may then solve for the only remaining perturbative held, h® x , in terms of 
J . Performing the same trick with equation (3.8) we see that the d x b®{x) disappears under integration, 
and the result gives us J = J{h® x ,E). Combining these two results we obtain an expression for the 
conductivity. To write the expressions compactly it is useful to introduce some notation; let 


h 


h 


h 


dx ( 6 — 87t T + 7T T d z log 



Qxx Qzz 

Qtt Qyy 



(3.14) 


8 Here £ denotes the Lie derivative along the indicated field and i^ indicated contraction with the vector field. 

9 The derivation of the conserved charges, J, and, Q , is presented in full generality in the appendix of [12]. We have also 
explicitly checked that the derivation of the conserved currents and charges holds for the model we are considering. 
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In terms of these integrals we find that the DC electric conductivity and the thermoelectric conductivity 
are given as 


- 2/ i 

Q_ = 4 ttI 3 

TE 2 hh-ll' 


(3.15) 


The calculation of the thermal and thermoelectric conjugate conductivities, R, and a proceeds in a 
very similar fashion. In this case however our ansatz for the form of the perturbations contains two time- 
dependent components b x = b x (z,x) — et r A(z, x) and ht x = ht x (z,x ) — etQtt( 2 ,x) ; % 7 . 111 The restriction 
of regularity in ingoing coordinates again determines the leading scalings for the fields. The result for the 
htx and b x fields is found to be: 


h tx (z, x) ~ . (h° tx {x) + 4nTh l tx (x)(l 
V "w K 

b x (z, x) ~ b x (x) + 0(z - 1) 


z) log (47 tT(1 - z))^j + 0(1 


z) 


(3.16) 


while all other fields exhibit the same scalings and constraints as displayed in equations (3.6). The con¬ 
straints displayed in (3.7) are modified by the presence of a logarithmic term in the ht x expansion and the 
lack of an external electric field. The resulting constraints are given by: 

h °z X (x) = htx(x), h° zz (x) = 2h 0 tz (x)-h 0 tz (x), b° z (x) = b° t (x), h\ x = ( 3 - 17 ) 

In order to obtain expressions for the response functions which do not involve the background fields 
we expand the conserved quantities of, J, and Q to leading and subleading order, respectively, near the 
black hole horizon. The leading order term in the expansion of Q again indicates that h^ x (x) is: 

Q = — 47t T h® x (x) (3.18) 


While the remaining two constraints displayed in equation (3.19) may be integrated over a single period 
such that they become linear equations in Jh^ x and r. 



A 


h tx \ 

°Qzz) 


+ b° t '(x) = 0 



h° tx 

2irTQ zz 


7tT< 9 2 log 


QxxQzz 

QttQyy 


871T + 6 


JA 0 \ 
4tt TQ zz ) 


+ 


d hj Z 
dx Qzz 


+ r = 0 


(3.19) 


10 The form of the coupled metric and gauge field perturbation is chosen such that the equations of motion remain time- 
independent. In writing the form of the metric perturbation as above we have tacitly modified our background metric ansatz 
such that Q xz = 0. This is done to avoid the more complicated form of the metric perturbation necessary if Q xz ^ 0. It 
has been checked that the response function results derived from the near horizon behaviour do not change as a result of this 
modification. 
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This information is sufficient to calculate the thermal and conjugate thermo-electric conductivities: 


_ Q _ 167 t 2 T I 2 

' = 2 hh-il 


J _ 4 tt/ 3 

Tt 2 hh-il 


(3.20) 


Once we have an expression for k we may easily calculate the conjugate thermal conductivity, n which 
corresponds to the the thermal conductivity at zero electric current. This is done via the relation 


k = k — 



8tt 2 T 

“tT 


(3.21) 


It is reassuring to note that the a = a as it should since our system does not break time reversal 
invariance. We may also check that the high temperature limit of the response functions behaves in the 
expected fashion. This limit may be extracted by expanding around — = T 0 _ —> oo. In this limit the 
solution resembles that of a Schwarzschild-AdS black hole and so appropriate field substitutions are: 


Ao —> S e fi. Q zz —> 1 + S e A Q zz , Qtt —> 1 + <5 e AQtt, Qxx —)• 1 + A Q xx , (3.22) 

Qyy ^ 1 T ^Qyy Qxz ^ A Qzz: 4* ^ A(^>, M ^ b 

where we may now expand around 5 e = 0. 11 We find in this limit that a = a —> 0 and R —> oo diverges 
as the effects of momentum dissipation are removed from the system. The electrical conductivity, a, 
asymptotes to Z( 0) = 1, the known result for Schwarzschild-AdS 4 . Furthermore, we can confirm that in 
the hydrodynamic regime the relation between the conductivities (3.4) is satisfied. In fact, from (3.15) and 
(3.20) we learn that in the hydrodynamic limit the background geometry should satisfy 2/i I 2 — if = 0. 1-> 

It was shown in [28] and [12] that the breaking of translational invariance in holographic models 
produces a low frequency AC conductivity well fit by the Drude model of conductivity. We confirm that 
this is true in our model in the appendix. It will therefore be interesting to test if our model obeys the 
Wiedemann-Franz law at any point in its phase diagram. We therefore calculate the Lorenz factors for 
both k and R as follows: 

f _ 87r 2 I 2 

h 

T 4 tt 2 (2 hh-il) 

n 

In §4 we will see that the that these expressions are not constant in either the metallic or insulating phases 
of our model, confirming departures from the Wiedemann-Franz law. 

11 The high temperature limit may also be accessed by undoing the scalings of the action and equations of motion which 
are used to keep the location of the horizon fixed at z = 1 as the temperature is changed. As the location of the black hole 
horizon, z p , then approaches the conformal boundary as ^ —> 00 we may find the high temperature behaviour of the response 
functions by expanding around z p = 0. The appropriate expansions of the fields in this limit is determined by their known 
conformal behaviour. As expected the results from this approach agree with those presented in the text. 

12 Strictly speaking this relation is valid in the translationally invariant case, where the DC conductivity contribution of 
the Drude peak diverges. 


(3.23a) 

(3.23b) 
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4 Transport results for holographic systems 

We regard the present effort as a first step of exploring the vast phase diagram of the the effective holo¬ 
graphic theories, identifying interesting corners for further study. In this section we discuss features of the 
DC conductivities and the optical conductivity, and qualitative changes in the physics as we vary param¬ 
eters. We argue that these changes indicate the existence of metal-insulator quantum phase transitions at 
various loci in the phase diagram. 13 

Let us start with a brief description of the numerical techniques we employ, before describing our main 
results. For those interested, much more detail regarding the development and testing of our numerical 
methods is provided in Appendix A. In order to solve our PDEs numerically we discretize them using 
spectral methods [7, 33]. A Chebyshev grid was employed in the radial direction and a Fourier grid in 
the spatial direction, which imposes spatial periodicity, as discussed above. The solution of the non¬ 
linear background equations uses both Newton and quasi-Newton methods, whereas the solution of the 
perturbation equation only necessitates the inversion of a matrix, for which we use direct methods. 

4.1 DC conductivities 

We start by exploring the direct electric and thermoelectric conductivities. These are calculated using 
horizon data, as explained in §3. The electric conductivity is also calculated as the zero frequency limit 
of the optical conductivity. Besides providing a check of the numerics, the low frequency behaviour of the 
optical conductivity helps in understanding and elucidating the IR physics. 

Since we are working at finite temperature, we cannot probe the metal-insulator quantum phase 
transition directly. Nevertheless, by interpolation of our knowledge of both the temperature dependence 
of the DC and optical conductivities to sufficiently low temperatures, we can diagnose the presence of such 
transition as function of parameters. In our exploration we fix C = 1.5, and discuss the phase diagram 
as function of v as well as the perturbation wavenumber k. Varying these two parameters we find both 
metallic and insulating regimes, and transitions between them. 

In Fig. 1 we show a representative sample of the DC conductivities with k = 1 and varying v. We 
clearly observe the transition from an insulating to metallic behaviour as the value of v is increased. For low 
value of v, which includes the Einstein-Maxwell model of [28], we find a distinct insulating behaviour: the 
conductivity decreases at low temperatures and seems to vanish at zero temperature. On the other hand, 
for sufficiently large v the opposite behaviour is manifest: the conductivity is increasing with temperature 
and seems to diverge at zero temperature. In the transition region, the conductivity shows no distinct 
trend - this is the region which is a bad or incoherent metal at low temperatures. 

A similar trend can be seen in the DC thermoelectric conductivity, which is also monotonically in¬ 
creasing as we lower the temperature, for small values of v. As we increase v the curve begins to kink 
downwards at low temperatures until eventually a well defined turning point is formed. This turning point 
migrates towards larger temperatures as we continue to increase our control parameter. We conclude that 
the transition between metallic and insulating behaviour exists in this observable as well. 

13 The phrase, phase diagram, here refers to changing both the sources in a given theory (by tuning the period of the lattice 
set by k) as well as explorations across theories (by changing Lagrangian parameter v). 
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DC conductivity as a function of T 



DC thermoelectric conductivity as a function of T 


T/H 



Fig. 1: Plots of the DC electrical and thermoelectric conductivities for a range of models against T for various theories 
parameterized by v, with C = 1.5, k = 1 held fixed. We clearly see the existence of the metallic and insulating 
regimes separated by an intermediate region. 


It is interesting to note that no such transitions occur in the thermal conductivity, which displays a 
simple monotonic increase as a function of temperature for all values of parameters we examined. Thus 
our theories are all good thermal conductors. In some sense this is not surprising as the Wiedemann-Franz 
law is explicitly violated in our expressions for Lorenz factors in equations (3.23a) and (3.23b). This is 
confirmed by the numeric results presented in Fig. 3. 

4.2 Metal-insulator transitions 

We have seen that the finite temperature results lend themselves to interpretation as indicative of a zero 
temperature metal-insulator transition. We characterize such a phase transition by a qualitative change in 
the low temperature behaviour of the response functions. One can delineate regions of the phase diagram 
where the holographic theories we consider describes good metals, bad metals and insulators. We now 
search for such transitions as function of v and k. 

The result, displayed in Fig. 4, separates the parameter space into the regimes of clear metallic and 
insulator phases, separated by intermediate regimes in which the conductivity is neither monotonically 
increasing or decreasing. We see that increasing v reduces the relevance of the sourced inhomogeneity, 
such that the transition to a metallic phase occurs at lower values of k. 

The intermediate phase which straddles the phase transition region can be viewed as one where the 
competition between metallic and insulating orders is strong. It is interesting to speculate in analogy 
with the domain model for magnetic phase transitions, that one is encountering pockets of the competing 
phases. This leads to incoherence in the transport, reducing for instance the conductivity from its metallic 
value. 

We note that the phase separation between metallic and insulating behaviour occurs along a locus 
k+(v) which is monotone decreasing - as we increase the lattice wavelength, we encounter a transition at 
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DC thermal conductivity as a function of T 



DC thermal conductivity at 
zero electrical current as a function of T 



Fig. 2: Plots of both forms of the thermal conductivity, k and R. We note both the qualitative similarity of the 
two quantities and the insensitivity to the variation in the v parameter. In all cases that we have examined, 
including those associated with the data used to construct Fig. 4, the thermal conductivities were seen to increase 
monotonically with temperature. 


Lorenz factor as a function of T 



Lorenz factor at 

zero electrical current as a function of T 



Fig. 3: Plots of the Lorenz factors associated with k and R. The fact that these factors are neither constant as a function 
of temperature nor v indicates that the Wiedemann-Franz law is violated. This is in accord with the lack of a 
phase transition in the thermal conductivity. 


lower values of v. This suggests that the efficacy in translating the lattice between UV and IR regions 
in the geometry is playing a role in the presence/absence of “charge carriers”. 11 The DC conductivity 

14 We use the phrase “charge carriers” somewhat loosely since we are talking about transport in a strongly interacting 
system with no obvious quasiparticles. 
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Metal-insulator transition of DC conductivtiy 


1-2 o ooooooooooooooo 

ox XOOOOOOOOOOOOO 
>: X x x XOOOOOOOOOOO 
>c x x x x x xOOOOOOOOO 
111- + + + + + + x ooooooo 

-- + + + + + + + + + XOOOOO 

-- + + + + + + + + + + OOOOO 

-- + + + + + + + + + +XOOOO 

.-- + + + + + + + + + + + XOOO 

l: - + + + + + + + + + + + x x O O 

-- + + + + + + + + + + + + XXX 

-- + + + + + + + + T + + + + + X 


o o o o 
o o o o 
o o o o 
o o o o 
o o o o 
o o o o 
o o o o 
o o o o 
o o o o 
o o o o 
o o o o 
o o o o 
x x o O 


0.9 i; 


+ + + + + + + 

+ + + + + + + 

+ + + + + + + 


+ 

+ 

+ 


0.8 i; 


0.7 i; 


+ + 
+ + 


+ + 
+ + 


0.01 0.02 0.03 0.04 0.05 0.06 0.07 

V 


0.08 0.09 0.1 


Fig. 4: The phase plot in the (■ v,k ) plane illustrating the metal insulator transition. The green data points correspond 
to regimes which are clearly metallic in character with monotonically increasing DC electrical conductivity at low 
temperatures. Likewise the red data points correspond to insulating phases characterized by a monotonically 
decreasing DC electrical conductivity at low temperatures. The intermediate region between these corresponds 
to the transition region where a turning point is still evident in the profile at low temperatures. 


is effectively a proxy for the weight of the charge carrier spectral function at vanishing frequency. The 
support of this spectral function is localized in the vicinity of the horizon. This immediately follows from 
the fact that we have a membrane paradigmesque formula for the conductivity. 

Finally, we specialize to v = 0 which describes the Einstein-Maxwell model of [28]. This allows us to 
probe the existence of a phase transition as a function of the wavenumber k alone. 15 The results for the 
DC electrical conductivity are displayed in Fig. 6 and clearly show the existence of a quantitative change 
in the temperature dependence as k is varied. However at temperatures which we may reliable access the 
phase is always insulating for our choice of parameters (with C = 1.5 held fixed). This is in contrast to the 
metallic phases observed for the parameters chosen in [28]. 16 We have seen evidence that the transition 
may occur at lower temperatures however these are difficult to probe reliably. The qualitative shift in 
behaviour when moving from finite v to v = 0 can be illustrated by direct examination of the background 
field solutions as seen in Fig. 5. 

4.3 Optical conductivity 

It is interesting to examine the AC electrical conductivity in the vicinity of the transition between good 
metals and insulators. We observe that the influence of the v parameter on the profile of the real and 
imaginary parts of the AC conductivity is minimal until — is sufficiently small. This is in keeping with 
the profiles of the DC electric conductivity shown in Fig. 1, where the profiles do not begin to strongly 
differentiate until values of — ~ 0.06 are reached. At lower values of — the characteristic profiles associated 

/X /X 

15 Metal-Insulator transitions in Einstein-Maxwell theory deformed by helical lattices were discussed in [15]. 

16 We have checked that for the parameters C, k chosen in [28] we also encounter metallic phases. 
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Gauge Potential 


T=0.0015, k=2.02, u=0.002 


Neutral Scalar 



Real Qxz 


0 

-i 



T-0.001 5, k-2.02, u_0 Neutral Scalar 



Real Qxz 



Fig. 5 : Setting v to be strictly zero means that the IR evolution of the scalar field flattens out. As this behaviour of the 
scalar controls many aspects of the IR physics qualitative changes in the response functions are to be expected. 

with Drude behaviour begin to acquire a distinct spread as a function of v. This AC conductivity of 
manifestation of the metal-insulator transition, cf., Fig. 7. 

Our testing the zero frequency limit of the optical conductivity, agrees cleanly with the explicit evalu- 
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Metal-Insulator phase as a function of k 



Fig. 6: The DC conductivity in the limit v = 0 for different lattices. We notice that for lower values of the wavelength 
k the conductivity is decreasing at low temperatures. As we increase k the DC conductivity begins to flattens at 
lower temperatures though for this choice of amplitude, C = 1.5, for the sourced inhomogeneity the conductivity 
never transitions into the metallic phase for the temperatures which we can reliably access. 


ation of the DC conductivity using the membrane paradigm formulae. In the real part of <r(tu) we observe 
the characteristic Drude peak. We can confirm that the intercept cr(to = 0) agrees with the DC result 
obtained from the membrane paradigm, providing us with a nice consistency check of the numerics. 

Inspection of the low frequency optical conductivity lends further evidence to our physical picture of 
the low temperature transport. In Fig. 8 we confirm that as the temperature is lowered for parameter 
choice in the insulating phase spectral weight is shifted from lower to higher frequencies, as expected. 

We also check for potential mid-IR scaling regimes. A convenient way to identify this behaviour is 
using the diagnostic quantity which we label F(tt>) (cf., [14]): 


F{ to) 


1 + tt> 


o~o)r 

k(»r 


(4.1) 


If |<r(tr>)| develops a scaling regime such that it behaves as |cr(tt>)| ~ C\ + C 2 tv 1 ", the quantity F( to) will be 
equal to a constant given by the scaling exponent v. 

In Figs. 9 and 10 we plot F(tt>) for a variety of values of ^ and v, including those displayed in Fig. 7. 
We observe that as j l is decreased the profiles of the diagnostics F(\v) begin to flatten out and the existence 
of a scaling regime becomes a real possibility. However the existence of such a scaling regime requires the 
fine tuning of v (or 1) to a very narrow parameter range. This ranges questions regarding the generic 
nature and therefore physical significance of the behaviour. However, even when scaling regimes do exist, 
the scaling exponent seem to be generically different than 2/3. Further questions regarding the significance 
and robustness of this scaling regime is postponed to future work. 1 ' 

17 For a discussion of the difficulties involved in working in this frequency regime please see Appendix A. 
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Re(o) Re(o) Re(a) 


Real conductivity at T|i=0.015949 



Real conductivity at T4i=0.012167 



Real conductivity at T/p=0.0085301 



Imaginary conductivity at T4i=0.015949 



Imaginary conductivity at T4i=0.012167 



Imaginary conductivity at T/n=0.0085301 



Fig. 7: We examine the effect of the control parameter v on the AC conductivity as ^ is lowered. As the value of ^ is 
lowered the influence of v on the form of the real conductivity curves becomes more apparent. In this limit the 
curves for different values of v can be seen to clearly differentiate at lower values of in. 

4.4 High Temperature Limit 

As described in §3 the form of the response functions should exhibit high temperature behaviour consistent 
with that of the Schwarzschild-AdS 4 black hole. In this limit a —> 1, a = a —» 0 and R —> oo. Therefore both 
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Re(a) Re(a) Re(a) 


Real conductivity at v=0.01 


Real conductivity at v=0.03 




Real conductivity at v=0.05 Real conductivity at v=0.06 




Real conductivity at v=0.07 Real conductivity at v=0.1 




Fig. 8: Examining the real part of the AC conductivity as a function of ^ for a variety of v. At small v the constant 
curves intercept and overlap for lower values of ro. This can be interpreted as the AC representation of the 
incoherent phase. As v is increased the constant — curves separate and differentiate as we transition into a 

f 1 

conducting phase. 


the metallic and insulating low temperature phases must transition to that generic behaviour, determined 
by the conformal invariance of the UV theory, as the temperature is increased. 
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Mid-IR scaling test Mid-IR scaling test 


Mid IR scaling test at T/|i=0.019889 


Mid IR scaling test at T4i=0.015949 



Fig. 9: The behaviour of the diagnostic function F(tt>) as we scan for mid-range scaling behaviour as a function of v. 

At low temperatures, and for appropriately chosen values of v, the existence of a scaling regime is possible. For 
example, at temperature of — = 0.01267 we observe that the our test function flattens out at v ~ 0.09. 


An interesting point to note is that, as seen in Fig. 1, the DC electrical conductivity increases to values 
well above unity in the insulating phase. This can be understood as consequence of the sum rule obeyed 
by the optical conductivity - the suppressed low frequency spectral weight in the insulating phase has to 
transfer to high frequencies, which are still low compared to the scale of the chemical potential. Therefore 
for those insulating phases, a turning point must exist when the sign of the slope of the DC electric 
conductivity reverses, as the low temperature physics begins to transition to the generic high temperature 
behaviour. Conversely, the existence of such turning point could be another indication of the transition to 
an insulating phase. 

In Fig. 11 we plot the values of T cr i t = — t at the turning point versus v and k. For the DC electrical 
conductivity it can be seen that the position of this turning point decreases steadily as one moves from the 
insulating regime towards the metallic transition. We also plot the corresponding turning point in the DC 
thermoelectric conductivity (right panel). It should be noted however that a further, higher temperature, 
change of slope must exist for the thermoelectric conductivity, as it is expected to go to zero in the high 
temperature limit. 
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Mid-IR scaling test Mid-IR scaling test 


Mid IR scaling test at v=0.1 


Mid IR scaling test at v=0.09 



T/n=0.0085301 

T/p=0.012167 

T/(i=0.015949 

T/|i=0.019889 

T/|i=0.024003 




Mid IR scaling test at v=0.02 


Mid IR scaling test at v=0.01 



Fig. 10: The behaviour of the diagnostic function -F(id) as we scan for mid-range scaling behaviour as a function of 
temperature. The reversal of the slopes of F(m) as a function of temperature, for a fixed v, imply that it 
should be possible to tune to a scaling regime by appropriately specifying the temperature. 

5 Conclusions and outlook 

The main goal of this article has been a preliminary exploration of the large phase space of effective 
holographic theories, with a particular focus on understanding the efficiency of low energy transport in 
these models. More specifically we focused on the low frequency conductivity in a phenomenologically 
motivated AdS/CFT set-up and examined the predilection of the system towards metal-insulator phase 
transitions. We view our investigation as offering a large set of toy models, with features that can be 
dramatically different from those of the much studied Einstein-Maxwell theory. Our analysis has identified 
interesting loci in the phase diagram where qualitative changes of the IR physics take place. The natural 
next step is taking a closer look at these locations to get a better understanding of the physics that drives 
these transitions in holographic systems. 

One of the interesting regimes is that of incoherent, or bad metals which lie at the interface of the 
metal-insulator transition. A conjecture was put forward in [19], that transport in such incoherent metals 
was to be governed by diffusion processes. The models studied here provide a natural testing ground for 
this conjecture. 
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Turning point of the DC electric conductivity Turning point of the DC thermoelectric conductivity 




Fig. 11: The DC electric conductivity starts off at a finite value of r cr j t = — in the insulating phase and decreases 
rapidly as we approach the metal-insulator transition. On the other hand, in the thermoelectric conductivity, 
the existence of a turning point sets in as we enter the metallic phase and steadily increases in ^ as we tune 
v and k to move deeper into the metallic phase. 


Additionally, it is interesting to study the low lying quasinormal modes in the incoherent and insulating 
phases. The motion of those quasinormal modes in the complex plane can often elucidate the dominant 
physics governing the transition or cross-overs in the qualitative behaviour of the conductivity. In this con¬ 
text, it would also be interesting to quantify the scaling of the spectral weight and DC conductivity at low 
frequencies. In particular, examination of the residue of the lightest quasinormal mode (the hydrodynamic 
mode that contributes to charge diffusion), should illuminate whether the reduction in the conductivity as 
we enter the regime of incoherent metals (from the metallic side) is caused by the drop in spectral weight 
or if some other physics is responsible. 

From a gravitational viewpoint it would be useful to know if there is a characteristic feature of the 
near-horizon geometry which results in this phenomenon. Our preliminary investigations were inconclusive 
in ascertaining a sharp feature of the near-horizon geometry, which could be held responsible for the 
incoherence in charge transport. Ideally one would conjure up a geometric observable that is sensitive to 
transport. What we can definitely conclude is that for fixed v varying k demonstrates clear changes in the 
relevance of the inhomogeneity in all the metric components. Since decreasing k increases the relevance 
of the inhomogeneous scalar source, this makes it clear that smaller values of k will have a more stronger 
impact in the IR transport. Comparing across values of v is of course more complicated, since we exploring 
behaviour in the space of theories. Indeed even in the homogeneous case we see differing IR behaviour as 
we tune v. 

We have further seen that in certain small regimes of our parameter space intermediate frequency 
scaling may be possible. While our analysis of this effect has not been comprehensive, it appears that we 
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require a certain amount of fine-tuning in order to achieve scaling behaviour. Moreover, in most of the 
cases we looked at the exponent v was not 2/3, which was the value seen for the cuprates [34] and in v = 0 
models [27, 28]. While this is per se not surprising, since our explorations with r / 0 move us in the space 
of UV field theories, it is curious that the effect is rather dramatic in the low frequency conductivity. 

In conclusion, we hope that the results presented in this work illustrate that the class of models 
considered provides an interesting phenomenological environment for exploring detailed features of charge 
transport in holographic systems. We hope to continue this exploration, and address some of the interesting 
questions raised in the course of this work, in the not too distant future. 
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A Code Analysis 

In this appendix we provide further information on our numerical methods and exhibit some checks on 
the results presented in this paper. We first briefly describe our approach to solving the background and 
perturbation equations. We then present our convergence results and other checks of the numerics. 

A.l Numerical Procedure and Implementation Details 

The derivation of the equations, boundary conditions and gauge conditions were done in Mathematica. 
These expressions, in appropriately discretized form, were then exported to Matlab which served as the 
principal platform for equation solving and post-processing work. Significant portions of the code were 
transferred to C++ code, which utilized the Armadillo and BLAS linear algebra libraries. 

The non-linear nature of the background equations of motion requires an iterative process for the 
solution. This was accomplished using a combination of Newton method with line search and the quasi- 
Newton, Broyden method algorithms. The Broyden method was found to be the most efficient approach, 
as it did not necessitate the computationally expensive process of assembling the Jacobian matrix. As 
the results presented in this paper rely on knowledge of the temperature dependence of the conductivity 
at each point in the v, k parameter space, it was necessary to numerically solve the PDEs many times. 
As such the utilization of optimized code and efficient solver algorithms was important in maintaining a 
manageable computational load. An example of a useful solver strategy for the background was: 

• First attempt to use the Broyden method. If the initial guess is sufficiently close to the sought after 
solution, convergence should be reached rapidly. 
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• If the Broyden method fails to converge the solver switches to Newton method augmented with a 
three-point safeguarded parabolic line search. Once the norm of the residual has reduced below a 
safe tolerance, the Broyden method can once again be utilized to quickly bring about convergence to 
the required accuracy. 

• Once an initial solution is obtained suitably sized (adiabatic) variations in parameters allow for the 
mapping of parameter space, without having to resort to the Newton method. 

The solution of the perturbation equations is in some sense easier as the linear nature of the equations 
means that they can be inverted in one step without need for iteration. An additional complication is 
introduced by the fact that knowledge of the temperature dependence of the AC conductivity requires us 
to scan over both temperature and frequency at each value of the parameters v, k. As the quality of the 
numerical convergence is not uniform with temperature or frequency, we must be careful to establish which 
regions of the ( T,u,v,k ) space we can reliably probe and bear this in mind when analyzing our results. 
We also must be careful to ensure that we have sufficiently resolved the background solution relative to 
the desired resolution for the perturbation solution. 

A.2 Convergence Tests 

We now present the following convergence results for the background and perturbation solutions: 

• The background solution: 

— Convergence of the solutions to equations (2.2) as a function of the grid size N. 

— Convergence of the DeTurck gauge condition to zero as a function of N. 

• The linear perturbation solutions: 

— Convergence as function of the grid size N p for the linearized equations. 

— Convergence of the deDonder and Lorentz gauge conditions (2.11) to zero as a function of N p . 

— Convergence of the auxiliary (unimposed) horizon boundary conditions, as described below 
equation (2.13), to zero as function of N p . 

In Figs. 12 and 13, we consider the convergence of the solutions as a function of N, for various low 
temperatures and v. We consider the convergence properties of the solutions as a function of both the 
transverse grid size, N x , and radial grid size, N z , separately. This may be done by fixing one of the grid 
sizes and running convergence tests on the other. We display the results below for the convergence of the 
log of the norm of the difference in solutions for three separate, low temperatures. In the first set we fix 
N x = 45 and vary N z while in the second we do the reverse. From these tests we may draw the (perhaps 
expected) conclusion that convergence of the solutions depends more strongly on the radial grid. 

We note that the asymptotic expansion of the fields near the conformal boundary contain logarithmic 
terms at high orders of the expansion. Therefore the exponential convergence of the spectral methods we 
are using is expected to fail for fine enough grids. However, we find that for the range of parameters we 
consider here, this is not an important issue. Similarly, at sufficiently low temperatures we expect that 
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Log of norm of residues Log of norm of residues 


Convergence of norm of residues 


Convergence of norm of residues 




Fig. 12: Convergence of the log of the norm of the difference of the solutions as a function of N z for temperatures of 
(0.0016, 0.0021) for a variety of v with N x = 45. These temperatures correspond to ^ values of approximately 
(6.6e~ 4 , 8.6e -4 ) and therefore can be considered to be small on the scale set by the chemical potential. We 
see that we obtain exponential convergence as a function of N z , and that, while these plots become noisier for 
smaller temperatures, the exponential convergence reasserts itself as N z increases. 


Convergence of norm of residues 



Convergence of norm of residues 



Fig. 13: Convergence of the log of the norm of the difference of the solutions as a function of N x for the same values 
of v and temperature as the previous plot and with with N. = 45. We note that the convergence deviates 
significantly from exponential and that the small scale of the y axis indicates that convergence is occurring 
more slowly than for the case of increasing N z . 


finite differencing approximation near the horizon is more suitable. It may be that utilizing such methods 
together with domain decomposition approaches will allow us to reach still lower temperatures in future 
work. 

In Fig. 14 we perform a similar series of tests for the convergence of the norm of the gauge condition 
towards zero. We again find that better convergence behaviour is achieved by increasing N z in preference 
to N x . We also note that different scaling regimes may exist in the convergence of the norm of the gauge 
condition as a function of the resolution, cf., Fig. 15. 
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Log of norm of gauge Log of norm of gauge 


Convergence of norm of gauge 


Convergence of norm of gauge 




Fig. 14: Examining the convergence of the log of the norm of the gauge condition as a function of N z . For the lower 
temperature (leftmost graph) we observe exponential convergence with the now expected noise at lower values 
of We note that for the higher temperature (rightmost graph) larger magnitudes of v exhibit two distinct 
scaling regimes with a crossover occurring at approximately N z = 55. While the convergence is exponential in 
both cases it is markedly faster in one case. This may mean that it is necessary to go to higher resolutions if 
very accurate solutions are required in this region of parameter space. This issue however was not encountered 
in the results presented earlier in the paper. 


Convergence of norm of gauge Convergence of norm of gauge 




Fig. 15: For the same values of v and temperature as Fig. 14 we see that the convergence of the log of the norm of 
the gauge condition is significantly slower then exponential when N x is increased for a fixed N z . 


From this series of experiments we conclude that the best resolution for a fixed number of grid points is 
obtained when more points are allocated to the radial grid number, N z . As a verification of this approach 
in Fig. 16 we plot the log of the norm of the residues as above with N z and N x both increasing but with 
N z = N x + 20. 

We now consider convergence results for the perturbation equations. Firstly it was noted that asym¬ 
metric grids are less useful in this case and that the best results were obtained when N x and N z were 
increased in tandem. In addition it was found that for smooth convergence to be obtained the background 
should be available at a higher resolution then the perturbation resolution. We are now interested in the 
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Log of norm of residues 


Convergence of norm of residues 


Convergence of norm of residues 




Fig. 16: Convergence of the log of the norm of the residues with both N z and N x increasing but with N x lagging N z by 
20. Again the ^ values are approximately (6.6e -4 ,8.6e -4 ) moving from left to right. We see good exponential 
convergence even for the lower temperature case. 

convergence as a function of v, T and co. 

For our purposes there are two ro regimes where we must examine the convergence of the AC conduc¬ 
tivity. These are tt> <C 1, which is relevant for comparison to the DC conductivity, and 5 ~< tn <~ 10 
which is relevant for examining potential IR scaling regimes. In Fig. 17 we display some examples of the 
convergence of the norm of the difference of the perturbative solutions, and the norm of the gauge and 
auxiliary conditions as a function of N for tn = 0.1. In Fig. 18 we plot the analogous results for the real 
and imaginary parts of the conductivities themselves. We note that as yus lowered and v increased the 
convergence of the norm of the gauge and auxiliary conditions becomes more strained. We have however 
checked that for the data displayed in Fig. 1 both are on the order of 10~' 5 when the DC data points are 
read off. 

The situation is more complicated in the opposite limit which we examine in Figs. 19 and 20 where 
tn = 7. We note that while convergence is maintained, the rate of convergence for the auxiliary horizon 
constraints becomes problematic as the temperature is lowered. It is instructive to examine the form of 
these auxiliary constraints corresponding to the results displayed Figs. 9 and 10. This is done in in Figs. 22 
and 23. We see that the constraint violation worsens as the temperature is decreased or v increased, and 
that the violation is worst in the regime of 3 ~< tn <~ 5. It has been observed that as one adjusts 
parameters further into these regimes numerical artifacts appear in both the imaginary conductivity and 
the diagnostic function, F( tn). However, given that the conditions are under better control for tn > 5, we 
believe changed that the results displayed in Figs. 9 and 10 are qualitatively correct. This conclusion is 
supported by the observations that: 

• As seen in Figs. 19 all quantities, including the auxiliary horizon constraints, remain in a convergent 
regime down to temperatures significantly below the lowest used in the construction of Figs. 9 and 
10 . 

• Even at these lower temperatures the physical variables in which we are interested (the real and 
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Log of gauge condition norm Log of auxiliary condition norm Log of norm of the difference in solutions 


Norm of the difference at T4i=0.052109 to=0.0098 w/T=0.1 


Norm of the difference at T4i=0.0011499 co=0.00028 co/T=0.1 




Fig. 17: Semi-logarithmic plots of the norm of the difference in solutions, and the norm of horizon auxiliary conditions 
and of the gauge conditions versus N at a low and high temperature for ^ = 0.1. We note that while still 
convergent, the resolution of the auxiliary conditions must be monitored closely for small values of — and larger 
values of v. 

imaginary parts of the conductivity) exhibit sensitivity to N only up to a small percentage of their 
magnitudes. This can be seen in Fig. 20 by comparing the degree of fluctuation at larger values of 
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Imaginary Conductivity Rea l Conductivity 


Real conductivity at T/|i=0.052109 to=0.0098 w/T=0.1 


Real conductivity at T/|i=0.0011499 w=0.00028 co/T=0.1 



Fig. 18: Plots of the real and imaginary conductivities versus N for the same choice of parameters as in Fig. 17. We 
note that, as expected, convergence is slower for larger values of v. This is particularly evident for the imaginary 
part of the conductivity 

N to the scale of the y-axis. 

As a final illustration of numerical control in the mid-IR regime we plot, in Figs. 21, the real and 
imaginary conductivity, and auxiliary horizon constraints at a temperature slightly below the lowest tem¬ 
perature used in Figs. 9 and 10. We see that while the behaviour is qualitatively similar to the lower 
temperatures we have just described, the convergence rate of the auxiliary horizon constraints has im¬ 
proved substantially. In addition is is clear that the real and imaginary parts of the conductivity have 
effectively become independent of N by the time we reach resolutions comparable to those used in our sim¬ 
ulations (~ iV = 40). Future work in this direction may involve the use of alternative numerical methods 
(for example finite difference discretization) to tackle these numerically difficult regimes. 
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Log of gauge condition norm Log of aux j|j ary condition norm Log of norm of the difference in solutions 



Norm of auxiliary conditions at T/|i=0.027214 a>=0.406 (o/T=7 



Norm of auxiliary conditions at T4i=0.0025647 co=0.0434 co/T=7 




Fig. 19: Semi-logarithmic plots of the norm of the difference in solutions, and the norm of horizon auxiliary conditions 
and of the gauge conditions versus JV at a low and high temperature for = = 7. We note the unusual dual 
scaling regimes for the gauge conditions in the higher temperature case. More importantly we note the slow 
rate of convergence of the auxiliary conditions as the temperature is decreased. 
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Imaginary Conductivity Real Conductivity 


Real conductivity at T/|i=0.027214 ro=0.406 a)/T=7 


Real conductivity at T4i=0.0025647 co=0.0434 to/T=7 




Imaginary conductivity at T/|i=0.027214 co=0.406 w/T=7 

—v=o.ori 

—v=0.02 
v=0.03 
— v=0.04 
v=0.05 
v=0.06 
—v=0.07 
—v=0.08 
—v=0.09 
—v=0.1 



Imaginary conductivity at T/|i=0.0025647 w=0.0434 co/T=7 



Fig. 20: Plots of the real and imaginary conductivities versus N for the same choice of parameters as in Fig. 19. Again 
we note that increasing v makes the convergence more difficult. In this case the effect is most notable in the 
real part of the conductivity. 


A.3 Fit to Drude Form 


Another useful test, described in [28], is comparing the low frequency behaviour of the AC conductivity to 
the Drude form of the conductivity, expected on general grounds [22]: 

Kt 


cr(w) = 


( 1 . 1 ) 


1-iwr 

where K is a constant. This test is useful only in the metallic phase or above the critical temperature for 
the onset of the insulating phase. Therefore at moderate values of ^ the AC conductivity is well modelled 
by the Drude behaviour as demonstrated by the examples in Fig. 24. This may be interpreted as an 
important additional test of the numerics- the quality of the fit to the Drude form being indicative of good 
convergence of the perturbation equations over the range of uo examined. 

An additional important test of the Drude behaviour is that the coefficient, K, should match with the 
residue of the zero frequency pole in the imaginary conductivity obtained in the homogeneous case: 


K 

Im(o-homM) ->• —, 

(jj 


as cj 


0 


( 1 . 2 ) 
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Log of auxiliary condition norm Log of auxiliary condition norm 


Norm of auxiliary conditions at T4i=0.0076422 co=0.108 w/T=6 


Real conductivity at T/|i=0.0076422 co=0.126 co/T=7 





Fig. 21: An illustration that the auxiliary horizon constraints are well within a convergent regime at temperatures below 
those used in Figs. 9 and 10. It is also evident that the real and imaginary parts of the conductivity are well 
converged relative to their magnitudes at these temperatures. 

In order to perform this calculation we solved the ODEs which govern the homogeneous phases of this 
model at the same point in parameter space as our PDEs (2.2) and extracted the residue of the pole. 
These ODEs may be derived from the (gauged) PDEs via the following substitution: 

Qxz(z,x) ->• 0, Qzz(z,x)-> Pi(z), Qtt(z,x) -t P2(z), (1.3) 

Qyy(z,x) -t Q xx (z,x) -f p 3 (z) A 0 (z,x) —>• Aq(z), $(z,x)-+$(z) 

Once the ODEs have been derived the horizon conditions can again be obtained in the same manner 
described in §2 for the PDEs. In this case these boundary conditions relate the value of the fields and their 
first radial derivative at the horizon. The conformal boundary conditions can also be written as a simple 
mixture of Dirichlet and Neumann conditions as in the PDE case. 

Pi(0) = l, P2(0) = l, P3(0) = l, $'(0 )=A, A o (0) = p (1.4) 

The ODEs were solved via a simple adaption of the Chebychev grid spectral technique described in §4. 
The results also provided an additional sanity check on the solutions of PDE background equations (2.2). 
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Auxiliary condition norm Auxiliary condition norm 



Fit of auxiliary condition T/^i=0.019889 



Fit of auxiliary condition T4t=0.012167 



Fit of auxiliary condition T4i=0.015949 



Fit of auxiliary condition T/|i=0.0085301 


Fig. 22: The auxiliary horizon constraints corresponding to Fig 9. We see that while very well satisfied in the vo —> 0 
limit there is a regime of significant constraint violation prior to the return of greater accuracy in the mid IR. 

When these equations were solved with the sourced inhomogeneity turned off, $ 1 ( 2 :) = A, good agreement 
was obtained with the ODE solutions. 

Considering the linearized perturbation equations around these ODE backgrounds we find that it is 
sufficient to restrict ourselves to perturbations of the following form: 

h xt = h xt (z) e~ lult , b x = b x e~ lU}t (1.5) 

Utilizing the background equations of motion it can be checked that there are two independent equations of 
motion which can conveniently be decoupled from each other and written as: (i) a first order equation for 
h x t in terms of the background fields, and (ii) a second order equation for b x . We are only concerned with 
the equation for b x . Boundary conditions at the conformal boundary consist of b x (0) = 1 corresponding to 
a conveniently normalized external electric field perturbation in the dual QFT. At the horizon regularity in 
ingoing coordinates necessitates the following leading scaling of the field b x = V(z) b r x 9 . A suitable horizon 
boundary condition may then be obtained via expanding b r x g to leading order in (1 — z). 

Once the solutions of the ODE equations were obtained we extracted the coefficient of the pole and 
calculated a goodness of fit measure given by f 1 — J. The relevant results are displayed in percentage 
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Fit of auxiliary condition v=0.1 


Fit of auxiliary condition v=0.09 




-T/g=0.0085301 
T/g=0.012167 
T/|i=0.015949 
T/|i=0.019889 
-T/p=0.024003 



-T/H=0.0085301 
T/|x=0.012167 
T/|i=0.015949 
T/|i=0.019889 
-T/p=0.024003 


Fit of auxiliary condition v=0.01 


-T/(x=0.0085301 
T/g=0.012167 
T/p=0.015949 
T/|i=0.019889 
—T/n=0024003 



-T/H=0.0085301 
T/(i=0.012167 
T/|x=0.015949 
T/|i=0.019889 
-T/p=0.024003 


Fig. 23: The auxiliary horizon constraints corresponding to Fig 10. Again we observe the non-trivial dependence of the 
auxiliary constraints on the temperature, v and to with greater accuracy being obtained at lower and higher 
values of to. 


form in Table 1. A'i at is the result obtained for the inhomogeneous sources. 


V 

0.01 

0.02 

0.03 

0.04 

0.05 

0.06 

0.07 

0.09 

0.1 

—=0.052109 

0.57 

0.26 

1.08 

1.89 

2.69 

3.48 

4.26 

5.78 

6.53 

—=0.027214 

_ 

2.67 

3.53 

4.36 

5.16 

5.93 

6.67 

7.38 

8.7 

9.32 


Table 1: Testing the fit to the Drude form of the conductivity and checking agreement between the homogeneous and 
inhomogeneous solutions. 































Im(a) lm(o) 


Drude fit of imaginary conductivity as a function of temperature 


Drude fit of real conductivity as a function of temperature 




Drude fit of imaginary conductivity as a function of temperature Drude fit of real conductivity as a function of temperature 




Fig. 24: The fits of the Drude parameters t,K to the AC electric conductivities for two moderate values of ^ , 
corresponding to those listed in Table 1, for 0.01 > v < 0.1 . We see that the fits are good, particularly for 
lower values of the w where the convergence behaviour of the perturbation equations, and associated gauge 
and auxiliary conditions is best. In both cases the maximum ^ values plotted are ^ ~ 1.25. 
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